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In this paper we improve and extend a recent analysis which showed that the p-u mixing 
amplitude cannot be unambiguously extracted from the pion electromagnetic form-factor 
in a model independent way. In particular, we focus on the argument that the extraction 
is sensitive to the presence of any intrinsic ujj — > mi coupling. Our extended analysis 
confirms the original conclusion, with only minor, quantitative differences. The extracted 
mixing amplitude is shown to be sensitive to both the intrinsic coupling ujj — > tttt and to 
the value assumed for the mass of the p°-meson. 
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1. Introduction 

The p — uo mixing amplitude is traditionally extracted from the pion electromagnetic 
(EM) form-factor, F^q 2 ), as measured in e + e _ — > 7t + tt~ [|TJ. The non-perturbative strong 
interaction effects that produce the significant enhancement in the interaction around 
y/q 2 ~ 750MeV have been successfully parametrised using the vector meson dominance 
(VMD) model. Our interest lies in the appearance of the isoscalar uj meson in the isovector, 
p° meson dominated, pion form-factor. The traditional treatment neglects any intrinsic 
coupling of the u to the two pion final state (i.e., that not proceeding through a p° meson) 
and hence assumes that the p — u mixing amplitude is purely real. An argument due to 
Renard and others shows that, within certain approximations, the imaginary part 
of the p — uj mixing amplitude is cancelled by the intrinsic u decay contribution to the 
pion form-factor. This argument was critically questioned recently || and argued to be 
unjustified. In the following we discuss results from Ref. |3j], which extend and improve 
these arguments. We confirm the central conclusion, finding, in addition, a considerable 
sensitivity to the value assumed for the p mass. 
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2. The isospin pure basis 

Analysis of p— u mixing starts with the "isospin pure" fields, pi and Ui. These are by 
definition exact eigenstates of isospin. The photon to hadron interaction can be conve- 
niently described using a formalism in which the renormalised vector meson propagator 
is a 2 x 2 diagonal matrix. For the pion EM form-factor, in the pure isospin limit, we 
have no p-w mixing and no direct ujj — > tttt coupling and hence 

Fw = e^ Jpn f ^ ( (f £>L ) ( V ) ' ^ 

where e = |e|, / 7 y 7 is the coupling of the vector meson to the photon and g Pl7T7T is the 
coupling of the p to the two pion final state. Here and D 1 ^ are the scalar parts of 
the renormalised propagators for the isospin pure fields. As is standard in traditional 
treatments, coupling to conserved currents is assumed, allowing us to use D^q 2 ) = 
—g pu D(q 2 ), where D(q 2 ) is the "scalar propagator." 

The two eigenstates "mix" through the isospin violating mixing self-energy, I\. pul {q 2 ), to 
allow the decay uji — > pi — ► titt through off-diagonal terms in the dressed, isospin- violating 
vector meson propagator matrix ||. This introduces isospin violation (IV) in the vector 
meson propagator and we have (retaining only first order in isospin violation) EJ 



UJUJ 




+ 0((IV) 2 ). (2) 



However, a priori, in an effective Lagrangian model involving the vector mesons, all isospin 
violation sources are equally likely and there is no reason to exclude the possibility of the 
"intrinsic decay", uj — ► mr, through the coupling g Ul mz- The appropriate VMD-based 
expression for the pion form-factor in the isospin pure basis would then be given by 
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I h^l 1 n ($\ 

e q z — m z J {q z ) 

where a fourth term on the RHS involving both U pL0 and g^^n has been neglected since 
it is second order in isospin violation. One is always free to consider models where g W/7r7r 
is strictly zero, but there is no model-independent requirement that this be so. For 
the renormalised, isospin-pure propagators, D vv , we have used the physical p and uj 
propagators, since Dyy = D vv + O ((IV) 2 ) and since again we are working only to 
first order in isospin violation. For the physical p and uj propagators we will use here 
the usual Breit Wigner form with a momentum dependent width, i.e., Dyy = q2 _^i ^ , 

where m v (q 2 ) = m v —irhyTyhy{q 2 ) and where we have defined the momentum-dependent 
width Ty(q 2 ) = Tyhy(q 2 ) with hy(m v ) = 1 and where the mass parameter (my) and 
width parameter (IV) are real. 
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3. The physical basis 

The pion form-factor has two distinct resonance poles, yet Eq. @ has a rather compli- 
cated pole structure. Therefore, one might change to a basis that has only a sum of simple 
poles. This can be done by transforming to what is referred to as the "physical" basis. 
It has been traditional to assume that U puJ is a constant which allows the off-diagonal 
terms in the propagator to be removed by a rotation of the fields. However, recent work 
has argued that, when p—u mixing is generated only through vector mesons coupling to 
conserved currents, this amplitude is necessarily momentum dependent ||. This result 
is consistent with an earlier study which concluded that momentum independent mixing 
was acceptable for scalar particles, but not for vectors coupling to conserved currents || . 
Various demonstrative model calculations support this |/J. To make allowance for this, 
Maltman, O'Connell and Williams (MOW) defined the physical basis through two isospin 
violating parameters 0, p = pi — e\UJi, uj = oji + e 2 pi. This defines a transformation 
matrix, C, via 





= c = ' 7 1 Z ■ w 



u 

The off-diagonal element of the physical propagator was then defined in terms of the 
physical fields in the usual way through 

D^(q 2 ) = i J d 4 xe^(0\T(ff(x)^(0))\0). (5) 

As the vector mesons couple to conserved currents, we can replace D^ u (q 2 ) by —g^ u D{q 2 ). 
Eq. (|) leads us to 

DM) = DM) ~ tiDLia 2 ) + ^Jq 2 ) 

= ^ P (<7 2 )[rW<7 2 ) - eiP^ 2 ))- 1 + ^DUq 2 ))- l }DUq 2 ). (6) 

Note that D = CD 1 ^ . The physical basis is defined by requiring that there are no 
resonance poles in .D pu; (g 2 ), i.e., we choose E\ and e 2 such that D puJ (q 2 ) has no pole. The 
possible pole positions are the vector meson pole positions, m 2 and m 2 . These singularities 
can be removed if the numerator of D pu) (q 2 ) vanishes at these positions, leading us to, 

ei — — 5 o ' 62 ~ — 2 2" {') 

— nip — m p 

It should be noted that m 2 and w? w are the complex resonance pole positions. In the case 
where the mixing is momentum dependent we have e\ ^ e 2 such that the matrix C is 
not orthogonal and the transformation between bases is then not a simple rotation. In a 
similar fashion to Eq. (|5|), the coupling constants in the physical basis are defined via 

f — f _ f — \ ) 

J IP J IP I ^lj'yuji 9 pirn 9pii™- 

Note that ei^ r7r7r is second order in isospin violation and so is not retained. 
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In this way, MOW denned the form-factor in the physical basis by 
F 7r(q 2 ) = - [g^D wul f lul + g p7rn D pp f lp + g pnnDpujf -yu^} + background. (9) 

The non-resonant term, D pL0 , was then included in the background, along with any other 
non-pole terms of the Laurent series expansions of the propagators, to arrive at an ex- 
pression for the form-factor, 

F n (q 2 ) = Hie!, e 2 ) [P p + A(e u e 2 )e i ^ e ^Pj\ + background. (10) 

Here P p>u are simple poles and H is an overall constant. The quantities H, A and the 
Orsay phase, <ft, can be read off from Eq. (^j) as we will show later. In terms of the 
isospin-pure basis and the transformation matrix, C, Eq. (|9|) can be written, 



F n = -(f 7Pl f 1Wl )C T CD 1 C T C ( 9 n Pl 
6 y {Juji 




(11) 

Because of the closeness of the p and u pole positions, there is no practical distinction in 
their analysis between the two terms [i.e., Ii pu) {m 2 p ) and n pw (m^)], so we define a single 
parameter, e = e 2 = Il puj (m 2 p ) / (ml - m 2 p ) ~ ex, where Yl pLU (m 2 p ) ~ IL^m^,) ~ Xl^m 2 ) ~ 
11^(777^). In other words, from this point on we will assume (as is done in all standard 
treatments) that the momentum dependence of Tl puJ (q 2 ) is negligible in the vector meson 
resonance region. 



4. The Renard Argument 

The major conclusion of the MOW analysis concerns the competition between the 
two sources of isospin violation, namely, IL^ and g uinn . All isospin violation in the pion 
form-factor is usually attributed to p — uj mixing ||, because the intrinsic decay of the uj 
is assumed to be cancelled . Consider the ui pole term of Eq. (pf) . The coupling of the 
physical u to the two pion final state is given by, (7 W7rvr = g Ullxl , + eg Pl7T7T . It is certainly a 
reasonable approximation to assume that the two pion intermediate state saturates the 
imaginary part of U puj (q 2 ) around the p resonance region and that this is proportional to 
the two pion piece of the p self energy (which dominates the imaginary piece of the total 
p self energy, due to the strong p — > nn decay). We have then 

Im U pu) {m 2 p ) = Im W p Z{m 2 p ) = ^Im W p ;{m 2 p ) = -Gm p T p h p {m 2 p ) ~ -Gm p T p , (12) 

where we have defined G = 9uJI7r7T and assumed hJm 2 ) ~ hJm 2 ) = 1. We hence define 

Rpu{rn 2 p ) = ^(mj) - iGrhpTp. (13) 

Assuming saturation of the absorptive part by the two pion state, as described above, 
implies that 11^(777^) is real. Note that while the three pion state is kinematically acces- 
sible it also requires isospin violation and is, in addition, suppressed by the smaller phase 
space available. We follow standard practice and ignore it here. 
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Substituting Eq. fll3|) into the expression for e, we have 
n pw (m2) 



m 



2 _ 



iG- 



mi 



772 



2 _ 



-izT — zG, where z 



irhpTp 



mi 



rn 



2 _ 



T 



m- 



m p r p 



(14) 

~p ■■~u ---p ■•-w ""p " L p L p 

Note that since in the resonance region we are neglecting momentum dependence of the 
resonance widths, i.e. we take h^q 2 ) = h p (q 2 ) = 1, then we see that m^, m^, and z are 
all constants. Recall that here, since we are using just the single parameter e, the mixing 
amplitude Tl^m 2 ) should be understood to mean f[ pu} (q 2 ) with q 2 in the vector resonance 
region. Using this expression for e in the u pole term gives rise to Renard's cancellation 



9u 



9u 



mi 



G(l-z) + 



^-pU^ip) 



mi 



9 pi 



(15) 



when one makes the approximation z — 1. However, MOW find that z has a sizeable 
imaginary piece [[J z = 0.9324 + 0.3511 i, using the mass and width values from Bernicha 
e?; al. 0. For comparison, one finds z = 1.023 + 0.20387 using the values of Benayoun et 
al. Jl0| . The central observation of MOW was to point out that the deviation of z from 
unity leads to a substantial contribution to F^q 2 ) from g mw7r . 

At this point MOW neglected the e dependence of the constant, H(e\,e<z) in Eq. (]10[) 
and extracted G and n pw (m^) by fitting A and the Orsay phase 0. In order to do a 
more careful analysis of the G dependence of the extracted real part of the p — uj mixing 
amplitude, Yl puJ (m 2 ) we shall retain H(e) = ff(ei,e 2 ), which itself has a G dependence 
through e [via Eq. ([14])]. Let us write Eqs. (|9]) and (0), [or equivalently Eq. (§)] as 



-^7r(<? ) [f"fpi9pITT-IT-Pp fjUll^i.Pp Pu>)9 PITTK f-JUl P L>i7a>j7T7r] 



(16) 



where we have used the fact that e(P p — P w 



PpIL^P^. Defining the ratio 77 = 



/- 



we can use Eqs. ( |X4|) and (|TJ) to write Eq. fllED as 
Pr(<7 2 ) = -J^gp^APp + rjl-eiPp- P^ + GP^}} 

1 



7"J 

flpl 



-f^gpmn {P P (1 + ^r 7 f) + r 7 [-zzTP^ + zGP p + G(l - z)P,] } 



(17) 

Comparing Eq. (|TT| ) to our earlier Eq. ( JT0| ) it is seen to be a straightforward matter to 
identify the qualities H, A and the Orsay phase, <fi. This is the central result. 

Recall that to obtain this result, we followed standard treatments and neglected the 
momentum-dependence of the Breit-Wigner widths in the resonance region, h p (m 2 ) ~ 
hp(m 2 ) = 1. We also neglected the three pion loop contribution to p — u mixing, so that 
Im Ii pu {q 2 ) = Glm U pp (q 2 ) and assumed a constant value for Yl puJ in the vector meson 
resonance region. 

We arrive at a similar conclusion in a more straightforward way if we choose to stay in 
the isospin pure basis throughout and simply re-arrange Eq. @ to give 



Fn(q 2 



7P V 



+ 



m2 ,( ( l 2 ) Q 2 ~ m p(0' 2 ) 



ru(g 2 



9pi 



q 2 -ml(q 2 



+ 



im p T p h(q 2 



q 2 — m 2 (q 2 )' 



q 2 — rh? + irhpT p h(q 2 
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We see immediately that in the isospin pure treatment the G dependence (i.e., g W/7r7r 
dependence) is cancelled at q 2 = m 2 and somewhat suppressed around the pole region 
(where p — uo interference is most noticeable). Including the nonresonant off-diagonal 
propagator in the physical basis would make Eqs. ( |TTD and (|IBD identical, but as we can 
only fit resonant terms plus an unknown non-resonant background to data, there is no 
practical difference between these two equations. 



5. Results and conclusions 

decide upon Eq. (ITS), which we re- write as 



To determine values for II piJ j(m 2 p ) and G we must choose an appropriate form factor and 



—am 2 p 



q 2 — m 2 + irhpTp 



ft^ml) + G(q 2 - rh 2 p ) 

1 + Tj — 

q 2 -m 2 + im^T^ 



(19) 



We now perform a fit to pion form- factor data to extract values for m p , T p , Yl piJ j(m 2 ), G 
and the normalisation constant, a, using the fitting routine Minuit fll] (case A in Table 



0). We fix = 781.94 MeV and 1^ = 8.43 MeV, as given by the Particle Data Group 

It should be noted that we are here adopting this form to fit rather than using Eqs. (O) 
and (0), since we wish to isolate all of the G-dependence from the overall normalisation 
constant for the fit. It can be seen from Eq. ( |i~7|) that H in Eq. (|To| ) is G-dependent and 
so we can not use the fitting procedure adopted by MOW. By construction, the constant 
a in Eq. ([T9| ) is independent of G. For the analysis here we will assume the SU(3) value 
for 77 = g 7ti}I /g 7Pl , i.e., we will assume 77 = 1/3 |12| . However, from Eq. (|19D we see that 



all that can really be extracted from the analysis is riYl piAj and rjG and so if another value 
of rj is preferred our results can immediately be scaled in a straightforward way. Our 



data set consists of the 70 points listed in Ref. [13] in the region between 500 and 975 
MeV. While the preferred value for G is quite large - one usually expects isospin breaking 
at the few percent level, rather than 10% - the uncertainty is such that it lies only 2| 
standard deviations from zero. 

We can investigate the importance of direct isospin violation at the uo — > titt vertex by 
imposing the condition G = (case B). However, when fitting the p data, it is important to 
remember that the p parameters should not be process-dependent. As noted by Benayoun 
et al. since the p is a relatively broad resonance the value extracted for, say, the mass, 
can be greatly affected by the addition of other terms to a phenomenological form-factor 
[[Kj. To demonstrate that our conclusions concerning G and Yl puj do not depend on a 
particular choice of mass and width from one source, we have redone the fit using the 
Particle Data Group (PDG) values, which are averaged from a variety of processes |14 . 



In order to compare our results with the analysis of Maltman et al. || , we perform a 
fit to the same data using Eq. ( |T0D with H (e) treated as a constant (thereby absorbing it 
into the normalisation constant a). To first order in isospin violation one has 

Ae» = r(G(1 ~! ) ~^ f) . (20) 
1 + izrT + zrG 
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Parameter 


A 


B 


PDG 


MOW 


n pw (mp (MeV 2 ) 
G 
a 

m p (MeV) 
T p (MeV) 
X 2 /dof 


-6832 ± 1252 
0.102 ±0.04 
1.15 ±0.02 
764.1 ±0.7 
145.0 ± 1.7 
88/65 


-3844 ± 271 
(input) 
1.19 ±0.01 
763.7 ±0.7 
146.9 ± 1.4 
94/66 


-6298 ± 788 
0.065 ±0.033 
1.203 ±0.004 
769.1 (input) 
151.0 (input) 
172/67 


-6827 ± 1067 
0.102 ±0.04 
1.20 ±0.01 
764.1 ±0.7 
145.0 ± 1.6 
88/65 



Table 1 

Results from fitting our form-factor to data. 



The resulting fit parameters are shown in Table 1 under the heading MOW. We see that 
they are very close to the fit using the full expression Eq. flTTD, shown in column A. 

In conclusion, we have seen that, in agreement with the conclusions of Maltman et al, 
the pion form factor data supports equally well a large range of possible pairs of values 
for G and n pw . In other words, it is not possible to extract the p—u mixing amplitude in 
a model-independent way. The traditional method of extraction corresponds to assuming 
that there is no intrinsic u — > tttt coupling (i.e., that G — 0), which is highly unlikely. 
It should also be noted that these conclusions are entirely independent of what (if any) 
momentum-dependence is present in the p—uj mixing amplitude, since this was neglected 
in the resonance region in the usual way. 
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